Three-vortex configurations in trapped Bose-Einstein condensates 
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We report on the creation of three- vortex clusters in a 87 Rb Bose-Einstein condensate by oscillatory 
excitation of the condensate. This procedure can create vortices of both circulation, so that we 
are able to create several types of vortex clusters using the same mechanism. The three- vortex 
configurations are dominated by two types, namely, an equilateral- triangle arrangement and a linear 
arrangement. We interpret these most stable configurations respectively as three vortices with the 
same circulation, and as a vortex-antivortex-vortex cluster. The linear configurations are very likely 
the first experimental signatures of predicted stationary vortex clusters. 

PACS numbers: 03.75.Lm, 67.85. De 



Introduction — Quantized vorticity being a key feature 
of superfluidity, vortex dynamics and configurations have 
long been central topics in the study of Bose-Einstein 
condensates (BECs). In the context of trapped atomic 
BECs, vortex dynamics has enjoyed a resurgence of in- 
terest because of new physical effects associated with the 
trap geometry. Vortices have been produced in BECs 
by transferring angular momentum either by manipula- 
tion of the quantum phase and the hyperfine state of the 
atoms [1] or directly by mechanical means @, 0]. The 
majority of subsequent experiments have focused on ex- 
ternally rotated condensates, where vortices of the same 
circulation sign arrange into vortex lattices [4]. On the 
other hand, in BECs without external rotation, one can 
have vortices of both circulation signs coexisting. In a 
trap such combinations display fundamentally new ef- 
fects, such as stationary structures [1, 0, 0]. For exam- 
ple, in pancake-shaped BECs, a stable linear configura- 
tion exists for vortex tripoles if the interaction is large 
enough [y, [7| . This linear arrangement contains an anti- 
vortex at the trap center and two vortices equally spaced 
on either side of the anti- vortex. 

The purpose of this Letter is to report experimental 
evidence for the stability of specific three- vortex configu- 
rations in a trapped BEC. We create three- vortex states 
through an external oscillatory perturbation in a 87 Rb 
BEC held in a cigar-shaped magnetic trap. Our analy- 
sis of the three-vortex configurations clearly indicate the 
existence of stable linear structures. We interpret these 
as the vortex tripole structures of Refs. @, 0] described 
above. We also find structures where the three vortices 
form an approximately equilateral triangle. We interpret 
these as three vortices with the same circulation. 

In the absence of extensive experimental results on 
few-vortex dynamics and configurations in non-rotating 
trapped BECs, our results open up an intriguing direc- 
tion of vortex studies. By realizing the linear configura- 
tion in a cylindrical geometry, we demonstrate that the 



phenomenon of stationary vortex cluster configurations 
is robust well beyond the two-dimensional (2D) or pan- 
cake limit. In addition, the study of the dynamics and 
coexistence of vortices and anti-vortices is an essential 
ingredient for understanding evolution to turbulence as 
well as decay processes. More generally, interest in the 
interplay between vortex dynamics and confinement ef- 
fects extends beyond cold- atom physics, since analogous 
issues with magnetic vortex/anti- vortex effects in ferro- 
magnetic microstructures have attracted a great deal of 
attention in recent years [8]. 

Experimental results — The experimental setup and 
the vortex excitation procedure are described in detail 
elsewhere Briefly, we use a 87 Rb BEC with 1 x 10 5 

atoms held in a cigar shaped magnetic trap with frequen- 
cies ou z = 2tt x 23 and uj r = 2tt x 207 Hz. Vortices are gen- 
erated by applying on the BEC an oscillating quadrupo- 
lar magnetic field whose amplitude is much smaller than 
the trap fields. As a function of field amplitude we ob- 
serve an increasing number of vortices formed, up to the 
limit where a vortex tangle is observed [Tol ]. 

For this study, we focus on images where a clear 3- 
vortex configuration is observed within the condensate 
core. Therefore, we have restricted the amplitude and 
frequency of the excitation for generating only this kind 
of structure. The excitation time was fixed at 20 ms, fol- 
lowed by an equilibrating time of 20 ms before releasing 
the atoms for free expansion. This is important because 
different excitation times can lead to different physical 
situations. Vortices are produced while the BEC is still 
held in the trap. For our analysis we will assume that 
the in-trap structure is maintained during the 15 ms of 
time-of- flight. In fact, this is a widely used assumption 
in experimental vortex studies. 

We observe predominantly two types of configurations. 
The vortices are either distributed as a near-equilateral 
triangle, as in Fig. QJa), or as a near-linear array, as in 
Fig. [2(b). We analyze the geometrical distribution by 
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FIG. 1: Atomic optical density images showing configurations 
of stable vortices forming (a) an equilateral triangle, or (b) a 
linear array. Images were taken after 15 ms of free expansion. 



o 







> 

03 


0£ 



1.0 

0.9 

0.8 

0.7 

0.6 

0.5 

0.4 

0.3-I 

0.2 

0.1 

0.0 



52% 



60-100 100- 140 140- 



36.5% 



1 1 .5% 



60- 100 100- 140 140- 180 

a (degrees) 




FIG. 3: Observed relative frequency of 3- vortex configurations 
as a function of the angle a. The inset shows the expected 
distribution of a when the vortices are distributed at random 
positions. 



FIG. 2: (a) Definition of the largest internal angle a. (b) 
Schematics of the tripole configuration of vortices, arrows in- 
dicate the vortex circulation direction. 



recording the largest internal angle, a, of the triangle 
whose vertices are the vortex positions, as sketched in 
Fig. [2ja) . The histogram of Fig. [3] summarizes our results 
of more than 60 independent measurements, showing the 
relative frequencies of observed values of a. Two types 
of stable configurations have a clear prominence over the 
others, namely, an equilateral triangle type where the 
angle a is around 60°, and a linear array where a ~ 180°. 

Interpretation — In the inset to Figure 02 we show 
the expected distribution of the largest angle a when the 
three vortices are distributed at random positions within 
the visible BEC core, taken to be an ellipse with aspect 
ratio ~ 1.5. The intermediate a G [100°, 140°] arrange- 
ment would be expected to be almost as numerous as the 
equilateral-like a G [60°, 100°] arrangement. To explain 
the fact that we get much fewer intermediate configura- 
tions, we consider the dynamics of three vortices with 
the same circulation. For simplicity, our discussion and 
the simulations of Figure HI use a circular trap, but this 
is sufficient for the present discussion. 

Three vortices placed symmetrically around the trap 
center (a — 60°) in a nonrotating condensate simply 
precess around the center, maintaining the equilateral 
configuration. If the initial configuration is moderately 
distorted from this symmetric configuration, the dynam- 
ics of the triangle shape is also moderate, so that a does 



not vary too much. Thus, a significant fraction of the 
configurations starting with a G [60°, 100°] stay within 
this range during the precession dynamics (Fig. H]A.,B). 

However, when the initially created vortex positions 
are such that a G [100°, 140°], the triangle formed by 
their positions changes shape as the vortices precess, 
so that the oscillations of a takes it in and out of the 
a G [60°, 100°] range (e.g., Fig. HP). Thus the preponder- 
ance of near-equilateral configurations in our experiment 
reflect the fact that there are many a G [60°, 100°] con- 
figurations which are 'stable' in the sense of maintaining 
the triangle shape while precessing. 

While the equilateral arrangement is reminiscent of the 
Abrikosov lattice observed by many authors |ll|, a nd of 
the three- vortex case studied in Refs. [12] and [13], our 
measurements explore the rather different physics of non- 
rotating condensates. Our results reveal that in the sit- 
uation without external rotation the equilateral arrange- 
ment also plays a special role, through maintaining its 
shape during precession. 

From similar arguments (Fig. HP), if all vortices 
have the same circulation direction, we would expect 
a smaller frequency of near-linear configurations with 
a G [140°, 180°], compared to the random-distribution 
case of Figure [3] inset. In our experimental data, how- 
ever, the linear array is more frequent than the random- 
distribution prediction. This indicates that in the linear 
configurations the vortex circulations are not all of the 
same sign. 

Refs. [6, 7] have shown that a tripole configuration like 
the one shown in Fig. Efb) can be a stable configuration. 
For a pancake-shaped condensate, the vortex tripole is 
metastable when the interaction strength is above a cer- 
tain value 0. In this situation, the interaction strength 
is expressed by the effective 2D coupling parameter g, 
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FIG. 4: (Color online.) Evolution of the largest angle a, 
in Gross-Pitaevskii simulations starting from various three- 
vortex configurations in a circularly trapped two-dimensional 
BEC. Initial configurations are shown on right. 

defined as 



where a is the scattering length, N is the number of par- 
ticles, and a± is the oscillator length in the direction 
perpendicular to the plane considered to be quasi-2D. 
Following Ref. [7] the condition for the tripole config- 
uration to be stable is g > 108. In our experimental 
geometry, the perpendicular direction is one of the radial 
directions, a± = a r = ^/Ti/mou ri so that g ~ 200. Thus, 
from quasi-2D arguments, we expect the linear tripole to 
be metastable in our setup, which supports our interpre- 
tation that the predominance of linear configurations is 
due to the tripole of Refs. 0,0]. The fact that our setup 
is not really quasi-2D indicates that the results of Ref. [7] 
are valid well beyond the pancake-shape limit. 

It is expected [7j that the tripole survives for a time 
of the order of Uq 1 , which for our case is about 10 ms, 
well within the 20 ms of waiting time before releasing the 
atoms from the magnetic potential. Thus tripole decay 
has a chance to reduce the number of observations of 
tripole configurations, as discussed further below. In fact, 
for a very long waiting time after excitation, observation 
of collinear 3- vortex configuration seems to be scarcer but 
fluctuation does not allow us to be more quantitative in 
this regard. 

The experimental method used here to produce vor- 
tex states has basically the same probability to produce 
both circulation sign vortices. Such equal chances come 
from the fact that the oscillatory field introduced by the 
extra coils modulates the equilibrium position as well as 
the overall trap shape introducing features that result in 
a vortex state of both circulation signs, as discussed in 
Ref. [3]. The equal chances for vortex and anti- vortex 
formation implies that in many of the observed distri- 
butions a collection of two vortices and one anti- vortex 
can occur. Since such a combination has a metastable 



linear structure, this would explain the prominence of 
near-linear structures in Figure 03 

Relative abundance of configurations — If equal 
chances to produce vortices and anti-vortices are con- 
sidered, we would expect that the probability Pu ne to 
excite a mixture containing one vortex with certain sign 
of circulation and two vortices with the opposite one 
would be higher than the probability Ptriangie to pro- 
duce three vortices with the same sign. Considering 
the data presented in Fig. [3l we have Pu ne ~ 0.37 and 
Ptriangie ~ 0.52. From a simple probability consideration 
and assuming the production of a sequence of vortices 
as independent events, we would have the probability of 
having three vortices of the same sign as 0.25 while for 
having one of opposite sign among them as 0.75. That 
would produce a ratio 1 : 3 and not 1 : 0.7 as observed 
in our experiment. This mismatch is probably due to 
the complex dynamics occurring after vortex nucleation. 
If the surviving chances were equivalent, we would ex- 
pect more events showing tripole configuration than tri- 
angular configuration. Decay mechanisms may, however, 
modify completely such a expectation. Vortex tripole 
may have a vortex/anti- vortex pair annihilated. There- 
fore, this may be indicative of higher instability for the 
situation composed of two vortices and one anti-vortex. 
For times longer than ~ u^ 1 , the tripole tends to decay, 
for example through one vortex migrating to the borders 
of the condensate generating surface modes, and leaving 
behind a vortex dipole [7[ . 

Let us consider that the tripole configuration decays 
exponentially, with time constant uu^ 1 ~ 10 ms, to a 
dipole configuration. In that case, supposing that the 
tripole was generated, the probability of observing it af- 
ter the 20 ms of waiting time decreases by a factor of 
e -2 w 0.135. Let us additionally consider that the trian- 
gular configuration (with all vortices of the same sign of 
circulation) is preserved during the waiting time. Taking 
all this into account, if a configuration with three vor- 
tices is initially produced, the probability of observing 
the triangle remains as 0.25, the probability of observ- 
ing a tripole is 0.10 and the probability of observing a 
dipole is 0.65 (i.e., the expected configuration after the 
tripole decay). If from these observations we only con- 
sider the configurations with three vortices, we shall have 
that Ptriangie ~ 0.715 and Pune « 0.285. This is much 
closer to the observed values in our experiment. Since 
this is not the unique decay mechanism (vortex/anti- 
vortex pairing may also be important) the tripole con- 
figuration should be undercounted. That may explain 
the relative values in the histogram of Fig. OH 

Discussion — In summary, we have observed 3-vortex 
cluster configurations where two different types are the 
most frequent: the equilateral triangle geometry and the 
linear array. We have argued that the triangular clus- 
ter type is formed predomniantly by vortices with the 
same sign of circulation while the linear array is numer- 
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ous because of configurations where the middle vortex 
has opposite sign. 

The present experiments are loosely related to the work 
of Ref. [14[ where different vortex configurations were cre- 
ated by merging multiple BECs. Almost linear or trian- 
gle structures can also be seen in their data. Our results 
complement those of Ref. [14| , using a quite different vor- 
tex formation mechanism. 

Vortex-antivortex dynamics in confined geometries is 
of broader interest than superfluid gases, since similar dy- 
namical issues are currently under intense study in the 
magnetic microstructure context [8]. Within the field 
of cold atomic gases, our experiment features a distinct 
manner for creating vortices of both circulation in the 
same condensate. Our results highlight several poorly 
understood aspects of vortex physics, in particular in- 
volving non-equilibrium dynamics. Since the imaging is 
performed after the released condensates have expanded 
and have undergone aspect ratio inversion, there are sev- 
eral issues related to the expansion dynamics. It is not 
known in any detail whether or not the presence of vortex 
clusters affects the inversion dynamics, or how the vor- 
tex configurations are rearranged as the cloud undergoes 
anisotropic expansion. 

Another issue raised by the current experiments is 
the stability of the tripole configurations in a cylindrical 
rather than pancake-shaped condensate, since the the- 
oretical studies @, 0] concern flat condensates. The 
abundance of the linear tripole configuration seen here 
suggests that the scissor oscillations used to excite the 
condensate may stabilize the tripole configuration, per- 
haps counteracting the dynamical instability discussed in 
Ref. [7] . These effects are clearly of basic interest but are 
yet to be theoretically analyzed in detail. In addition, 
investigation of the dynamics of superfluids containing 
both vortices and anti- vortices may be a fruitful alterna- 
tive way to understand basic decay mechanisms of quan- 
tum turbulence and other poorly understood phenomena 
in atomic BECs. 
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